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Abstract 

Spectra taken by the Mars Global Surveyor Thermal Emission Spectrometer (TES) have been 
used to monitor the latitude, longitude, and seasonal dependence of water vapor for ove - one full 
Martian year (March 1999-March 2001). A maximum in water vapor abundance is observed at 
high latitudes during mid-summer in both hemispheres, reaching a maximum value of -100 
pr-pm in the north and -50 pr-pm in the south. Low water vapor abundance (<5 ?^ 
observed at middle and high latitudes in the fall and wmter of both hemispheres. There are large 
differences in the hemispheric (north versus south) and seasonal (perihelion vers « s aphelion) 
behavior of water vapor. The latitudinal and seasonal dependence of the decay of the northern 
summer water vapor maximum implies cross-equatorial transport of water to the 
hemisphere, while there is little or no corresponding transport during the decay of the southern 
hemisphere summer maximum. The latitude-longitude dependence of annually-averaged water 
vapor (corrected for topography) has a significant positive correlation with 

negative correlations with thermal inertia and surface pressure. Comparison of TES ^results w 
those retrieved from the Viking Orbiter Mars Atmospheric Water Detectors (MAWD) 
experiments [Jakosky and Farmer , 1982] shows some similar features, but also man J sl 8™ c ^ 
differences. The southern hemisphere maximum observed by TES was not observ * d * 
and the large latitudinal gradient in annually-averaged water vapor observed by MAWD does 

not appear in the TES results. 
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1. Introduction 

Although water vapor is a minor constituent in the 
atmosphere of Mars, it still plays an important role. 
Along with the annual cycles of CO 2 and dust, the 
annua! water cycle is one of the three main cycles 
that define the present Mars climate. Since the first 
detection of water vapor on Mars by Sptnrad ct al. 
[1963] there have been many additional observations, 
both from ground-based observers and from space- 
craft. Jakosky [1985] and Jakosky and ffaberle [1992] 
give good reviews of water vapor measurements. 

Since the time of the Viking mission, the defini- 
tive description of the seasonal and latitudinal depen- 
dence of water vapor has been that derived from the 
Viking Orbiter Mars Atmospheric Water Detectors 
(MAWD) [Former et a/., 1977; Jakosky and Farmer , 
1982]. That data set shows a global annual average of 
about 10 pr-/im of water with higher values over the 
northern hemisphere in the northern summer reach- 
ing a maximum of ~90 pr-/im near the north pole at 
about L t — 120°. Although interannual variability 
has been noted [e.g. Clancy ct al . , 1992; Sprague ct 
o/., 1997], most ground-based observations of water 
vapor have been more-or-less consistent with the gen- 
eral pattern observed by MAWD. Notable exceptions 
to the MAWD trend were the observation by Barker 
et ai [1970] of a maximum in water vapor abun- 
dance (reaching 30—40 pr-/im) in the southern sum- 
mer season compared to the very weak maximum in 
the MAWD data of 10-15 pr-^m at that season, and 
the observation by Barker of a somewhat higher wa- 
ter vapor column (15-20 pr-/im) at all seasons during 
the two Mars years 1972-1976 [Jakosky and Barker 
1984]. 

The Thermal Emission Spectrometer (TES) on- 
board the Mars Global Surveyor (MGS) has been in 
mapping orbit around Mars since March 1999, and 
now provides an excellent opportunity to examine the 
column abundance of water vapor over an entire Mars 
year in great detail. Previous results of atmospheric 
thermal structure and aerosol distribution from TES 
are given by Smith ct al. [2001], Pearl et al. [2001], 
Conrath ct ai [2000], Smith ct al. [2000a, c] and 
Christensen ct ai [1998]. In this paper we report 
the results of one Mars year of TES monitoring of 
water vapor. 

In section 2 we give a brief description of the TES 
instrument and the data taken during the mapping 
orbit. More detailed information is given by Chris- 
tensen et ai [1992; 2001]. In section 3 we describe 


the algorithm used for the retrieval of water vapor col- 
umn abundance. In section 4 we present our results 
and compare them with Viking MAWD and ground- 
based observations. We discuss our findings in section 
5, and summarize in section 6. 

2. Data Set 

The Mars Global Surveyor mapping orbit gives one 
narrow strip of observations running roughly north- 
south. One day of data gives two sets of twelve such 
strips spaced roughly 30° apart in longitude, with one 
set taken near a local time of 0200 and the other near 
1400. For this work, we use only the nadir-geometry 
daytime (1400 hours) data. We do not use night- 
time data because they have lower signal-to-noiso and 
the atmospheric temperature profiles retrieved from 
nighttime TES data potentially have relatively large 
uncertainties near the surface where water vapor is of- 
ten concentrated. This is because the km-scale verti- 
cal structure of the near-surface boundary layer night- 
time temperature inversion observed in the MGS ra- 
dio occultation profiles [Hinson ct ai , 1999] cannot 
be resolved by TES. 

To obtain a global- and seasonal -scale view of wa- 
ter vapor abundance we average the long wavelength 
portion of daytime nadir-geometry TES spectra into 
bins that are 2° wide in latitude and 2° wide in sea- 
son ( L t ). All longitudes are binned together, so these 
represent zonal averages. This results in --1 1000 bins. 
Spectra are averaged separately for the two spectral 
resolutions (6.25 and 12.5 cm” 1 ) that are used by 
TES. The results presented below were retrieved us- 
ing data starting with the beginning of the mapping 
portion of the MGS mission at L s = 104° (March 1, 
1999) and extending for just over one full Mars year 
to i 3 =. 137° (March 28, 2001). The TES operated 
nearly continually during this entire period with only 
minor exceptions. Several million spectra were used 
in this analysis. 

3. Water Vapor Retrieval Algorithm 

Figure 1 shows a line-by-line synthetic water va- 
por spectrum convolved to the TES resolution of 6.25 
cm At this spectral resolution, TES has a sam- 
pling interval of about 5.3 cm -1 . To retrieve water 
vapor abundance, we use the rotation bands in the re- 
gion 240-360 cm” 1 (28-42 /im). At 6.25 cm -1 resolu- 
tion the hundreds of individual water vapor lines that 
be in this spectral region are combined into five main 
bands spaced roughly 20-30 cm” 1 apart. This spec- 
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tral region is a favorable place to perform a retrieval 
because at these wavelengths the TES signal-to-noise 
is relatively high (-500), the opacity from dust and 
water-ice is relatively low [Smith ct al. 2000a, b], and 
the distinctive, “picket-fence” spectral dependence of 
the water vapor bands allows the contribution from 
water vapor to the spectrum to be easily separated 
from the contributions from dust, water ice, and the 
surface. 

The TES contains a 2x3 array of detectors. Be- 
cause optical path lengths in the instrument are 
slightly different for each detector, the wavenumbers 
sampled by each detector for a given TES sample 
number are slightly offset (by as much as a couple 
of cm -1 ). Water vapor rotation bands have sufficient 
spectral structure (see Fig. 1) that spectra from dif- 
ferent detectors cannot be averaged together without 
significantly “smearing out” the bands. Therefore, 
we only use spectra from TES detector 2, which has 
the advantage of being one of the two detectors near- 
est the centra] ray of the instrument. We have cali- 
brated the wavenumber scale by shifting the grid of 
wavenumbers as a function of TES sample number 
until the shape and amplitude of water vapor lines 
observed in the TES spectra best match the synthetic 
spectrum. This wavenumber calibration is quite sen- 
sitive. A variation of one-tenth of the TES sampling 
(or about 0.5 cm -1 ) results in a noticeably worse fit. 
The resulting sampling of the water vapor spectrum is 
shown as points on the synthetic spectrum in Fig. 1. 

3.1. Model Atmosphere 

For each L , /latitude bin, the thermal structure of 
the atmosphere is taken to be the average of all at- 
mospheric temperature profiles available for that bin. 
Atmospheric temperatures are retrieved from each 
TES observation (each footprint of 2x3 detectors) as 
part of the standard operational TES processing us- 
ing a constrained linear inversion of observed radi- 
ance in the 1 5-/im CO 2 band [Conrath ct al. 2000]. 
The surface temperature for each bin is taken to be 
the average brightness temperature in the relatively 
transparent spectral window between 1285 and 1315 
cm -1 for all available spectra in that bin. 

Water vapor absorption is computed using the 
correla ted- k method [for example, Goodxj ct al., 1989]. 
Correlated-k coefficients for an array of different tem- 
peratures and pressures are computed using a line- 
by-line computation of water vapor transmittance. 
Water vapor line parameters are obtained from the 
GEISA-97 database [Hutson ct al. 1992]. For the 


water vapor lines of interest, the line parameters 
from GEISA-97 are essentially identical to those of 
the HITRAN-2000 database. An approximate correc- 
tion is made to the line widths to account for C 02 - 
broadening of water vapor (see section 3.5). The ver- 
tical distribution of water vapor is approximated as 
being well-mixed with the CO2 up to the water 
vapor condensation level and then zero above that. 
The condensation level is a function of both the at- 
mospheric temperature profile and the water vapor 
column abundance. Therefore, an initial guess for the 
water vapor column abundance is used to compute an 
initial condensation level and then the entire retrieval 
is iterated until a self-consistent solution is obtained. 
Typically no more than two or three iterations are 
required for convergence. 

3.2. Equivalent opacity and band depth index 

The first step in the water vapor retrieval algorithm 
is to convert the averaged observed radiance spec- 
trum, / 0 bsW, in each /latitude bin to an equivalent 

opacity spectrum, r 0 (v), through numerical solution 
of the radiative transfer equation: 

W«')=fW»Pl urf,*] c- T ^'» + 

/ flpXThvlfi-'^dr, (1) 

Jo 

where v is frequency, fi is the cosine of the emission 
angle, R[T, v) is the Planck function, T sur f is the sur- 
face temperature, T(r) is the atmospheric tempera- 
ture, and the integral is performed from the spacecraft 
(at r = 0) to the surface (at r = r 0 ). 

The above equation neglects scattering and the 
solar beam, and assumes a plane-parallel geometry. 
These are excellent approximations for the observa- 
tion of water vapor in the TES nadir- geometry mode 
at 240-360 cm -1 . The integral in opacity is evaluated 
by assuming that the opacity source is well-mixed up 
to the water vapor condensation level and then zero 
above that. 

The opacity spectrum obtained from solution of 
(1) contains opacity from water vapor, dust aerosol, 
water ice aerosol, and the effect of a non-unit sur- 
face emissivity. The contribution from water vapor 
is estimated by removing a continuum opacity from 
each of the five main water vapor bands. The un- 
derlying assumption is that the spectral signatures of 
dust, water ice, and surface emissivity are approxi- 
mately linear over the narrow spectral range (about 
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20 cm” 1 ) of each individual water vapor band [Smith 
et al 2000a,b; Band fie Id et al. 2000a]. An over- 
all water vapor band depth index for the 6.25 cm -1 
resolution data, rg 25 , is then formed by averaging the 
conti nu um- removed opacities for the five bands as fol- 
lows: 


r 6.25 = 


- r o(18)+^[7-o(19)+ro(20)+ro(21)]-r o (22)- 


To (24) + 2r 0 (26) - r 0 (28) - r 0 (29)+ 
2r o (30) - Tb(31) - r 0 (33) + 2r 0 (35)- 


r 0 (37) - r 0 (38) + 2r 0 (40) - r 0 (42) 


/ 12 , 


( 2 ) 


band depth indices. The data are well -described by 
the empirical linear relation 


7*6.25 = 1.45 ( 7 - 12.5 - 0.0035). (3) 

For each L $ /latitude bin, the 6.25 and 12.5 cm” 1 wa- 
ter vapor band depth indices are computed, and the 
band depth for the 12.5 cm” 1 resolution data is con- 
verted to an equivalent 6.25 cm” 1 resolution value 
using (3). The final band depth index for each bin 
is then given by the average of the band depth in- 
dex from the 6.25 cm -1 data and the equivalent band 
depth index inferred from the 12.5 cm” 1 weighted by 
the number of spectra at each resolution. 


where the numbers in parentheses indicate the TES 
sample number (see Fig. 1). Similarly, an overall wa- 
ter vapor band depth index can be formed for the 
12.5 cm” 1 resolution data, 7-12 5 , by averaging the 
continuum-removed opacity for the bands observed 
in 12.5 cm” 1 data. 

3.3. Cross-calibration of 6.25 and 12.5 cm” 1 
data 

Most of the spectra taken by TES before ^L s = 
90° (December 2000) were taken in the 12.5 cm -1 
spectral resolution mode, while most of the spectra 
taken after L, = 90° were taken in the 6.25 cm” 1 
spectral resolution mode. Figure 2 shows a compari- 
son of TES observations of water vapor at 6.25 cm” 1 
and 12.5 cm” 1 taken at nearly the same time (L s = 
307°) and location (79°S latitude). Although it is 
passible to perform the water vapor retrieval sepa- 
rately for each spectral resolution, we have found it 
more convenient to convert the strength of the wa- 
ter vapor bands observed in the 12.5 cm” 1 mode to 
an equivalent strength as observed in the 6.25 cm” 1 
mode, and then to perform a single retrieval using the 
equivalent 6.25 cm” 1 value of the water vapor band 
depth index. We have chosen to convert to 6.25 cm” 1 
resolution because the water vapor bands are much 
more clearly defined at 6.25 cm” 1 than they are at 
12.5 cm” 1 and so a more accurate retrieval can be 
performed. 

Most of the L t /latitude bins contain only one spec- 
tral resolution or the other, but there are still over 
1800 bins that contain at least 10 spectra from each 
resolution. These bins are used to find how the ob- 
served water vapor band depth at 12.5 cm” 1 relates 
to that observed at 6.25 cm” 1 . Figure 3 shows the 
relation between the 6.25 and 12.5 cm” 1 water vapor 


3.4. Water vapor column abundance from 
band depth index 

The conversion from a water vapor band depth in- 
dex to a column abundance is performed for each 
/> a /latitudc bin by computing a suite of synthetic 
water vapor spectra at a range of different column 
abundances. The synthetic spectra are computed at 
6.25 cm 1 spectra] resolution and at the wavenumbers 
sampled by TES. The synthetic spectrum is converted 
first to opacity spectra using (1) and then next to wa- 
ter vapor band depth index using (2). Water vapor 
column abundance is solved numerically by variation 
of the column abundance until the band depth index 
from the synthetic data matches the observed band 
depth. At this point a new water vapor condensation 
level is computed using the new estimate of water va- 
por column abundance and the entire retrieval process 
is iterated until convergence is obtained. 

Figure 4 shows the resulting fit for one L s /latitude 
bin (centered at 43° N latitude, L s = 85°) shown 
in terms of opacity. A continuum has been removed 
from the TES opacity spectrum. The fit is excellent 
over the five bands that are used in the retrieval (be- 
tween 240 and 360 cm” 1 ) and still remains good for 
the strong band centered at 225 cm” 1 and the weaker 
bands at 370-400 cm” 1 which are not used in the re- 
trieval. Both the relative amplitudes and spectra] de- 
pendence (or shape) of the bands are we II -reproduced. 
The retrieved water vapor column abundance for this 
case is 25 pr-pm. 

3.5. Water Vapor Line Broadening by 
Carbon Dioxide 

Synthetic water vapor spectra are computed us- 
ing the correlated- k formulation. The correlated-k 
coefficients are derived from a very high resolution 
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line-by-line computation of the water vapor spectrum 
using water vapor line strengths and widths from the 
GETS A- 97 database. However, the line widths given 
in the GEISA-97 database are for water vapor broad- 
ened by terrestrial air, not by CO 2 as appropriate for 
Mars. For water vapor abundances and atmospheric 
pressures typical for Mars, the net result at TES spec- 
tral resolution of an increase in line width is a decrease 
in water vapor transmittance (or an increase in water 
vapor opacity). It is important to carefully consider 
the effect of CO 2 broadening and to make this cor- 
rection to the line widths because the effect on the 
retrieved water vapor column abundance is substan- 
tial. The value for CO 2 broadening chosen below for 
this retrieval results in a reduction of ~40% in water 
vapor abundance. 

We will call the ratio of the CCVbroadening co- 
efficient for water vapor to the air-broadening coeffi- 
cient for water vapor, the “broadening ratio”. A value 
of the broadening ratio greater than unity indicates 
that water vapor lines are broader in a certain par- 
tial pressure of CO 2 than they are in the same partial 
pressure of air (as given by the GElSA-97 database). 
There is little experimental data on the broadening 
ratio. Laboratory measurements of CO 2 broadening 
in the near-infrared by Gamachc et ai [1995] indicate 
a broadening ratio between about 0.9 and 1 .8 depend- 
ing on the particular transition. However, none of 
the lines for which line widths measured were in the 
wavenumber range used in this retrieval. In modeling 
Venus near-infrared spectra, Pollack ct al. [1993] and 
Crisp ct al. [1991] have estimated CO 2 broadening of 
water vapor by using a constant broadening ratio of 
1.3 based on work by Howard ct al. [1956]. In model- 
ing Mars microwave spectra, Clancy ct al. [1992] have 
used a constant broadening ratio of 1 .6 based on work 
by Licbc and Dillon [1969]. None of these estimates 
for the broadening ratio is for the water vapor lines 
used in our retrieval. Given this lack of information 
and consensus, we choose to adopt a constant value 
for the broadening ratio of 1.5 that is roughly con- 
sistent with all of the above. We note (as elaborated 
below in section 3.6) that uncertainty in this param- 
eter leads to an uncertainty in the overall scaling of 
our water vapor column abundance estimates. 

3.6. Estimation of Uncertainties 

The algorithm used to retrieve water vapor column 
abundance has a number of steps, and there are at 
least seven main sources for uncertainty in the final 
retrieved value: 1) instrument noise and calibration, 


2) the retrieval of surface and atmospheric temper- 
atures, 3) the assumed vertical distributions of dust 
and water vapor, 4) the removal of a continuum opac- 
ity from sources other than water vapor, 5) the cross- 
calibration of 6.25 and 12.5 cm" 1 resolution data, 6) 
the calculation of synthetic water vapor opacity, and 
7) the estimation of the effect of CO 2 broadening on 
water vapor. 

The uncertainty introduced by instrument noise 
and calibration is relatively small. For daytime condi- 
tions, the typical instrument noise of 3x 10" 8 W cm" 2 
sr" 1 / cm -1 [Christensen ct al ., 1992, 1998, 2001; 
Christensen, 1 999] gives an opacity error of about 0.02 
for any one particular TES spectral sample. However, 
since we use 17 TES spectral samples (in the 6.25 
cm" 1 mode) the error in the observed band depth in- 
dex is much smaller. The absolute TES calibration is 
quite good, and errors that do occur tend to have lit- 
tle spectral structure (for example, a nearly constant 
offset). This type of error is effectively removed when 
we compute the band depth index and subtract the 
continuum opacity. Uncertainties contributed by the 
retrieved surface and atmospheric temperatures can 
be of order 0.03 for a particular TES spectral sample 
[Smith ct a/., 2000a], but this is a systematic error 
that tends to be nearly an offset in the opacity and 
is therefore removed in the continuum when the band 
depth index is computed. 

There is not sufficient information in a TES nadir- 
geometry spectrum to retrieve the vertical distribu- 
tion of dust or water vapor independently. TES limb- 
geometry spectra can potentially provide information 
on the vertical profile of dust and water vapor above 
a height of 10 km [Smith et a/., 2000c], but that re- 
trieval is not yet operational, and is beyond the scope 
of this work. We have taken the vertical distribu- 
tion of dust to be well-mixed, which has been found 
to be a good approximation for dust in Viking [Pol- 
lack et a/., 1977], Pathfinder [Smith ct ai, 1997] and 
TES data [Smith et ai , 2000a]. We have taken the 
water vapor vertical dependence to be that of a well- 
mixed profile below the water condensation level and 
zero above it. Analysis of the vertical distribution of 
water vapor by Rodin et ai [1997], Burgdorf et ai 
[2000], and Gunvell et ai [2000] are all consistent 
with this assumption. The water vapor profile could 
be made to follow the saturation condensation curve 
above the condensation level, but that value drops off 
very quickly with height and its contribution to the 
total column abundance is negligible. For the dust 
loadings observed by TES during the period consid- 
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ercd here, we estimate the uncertainty introduced by 
the assumed vertical distributions of dust and water 
vapor to be ~5% or less. 

The computation of a water vapor band depth in- 
dex by removal of a continuum opacity level contains 
the implicit assumption that the spectral dependence 
of the contribution to the effective opacity from ail 
other sources (dust and water ice aerosols and a non- 
unit surface emissivity) varies smoothly over the spec- 
tral width of one water vapor band (20-30 cm -1 ) in 
the spectral range considered (240-360 cm -1 ). The 
spectral dependence of dust and water ice aerosol do 
indeed have much less spectral structure on that scale 
than does water vapor [Smith et al., 2000a, b; Band- 
field et al., 2000a]. On a global scale, the two ma- 
jor surface emissivity spectral shapes (a basaltic and 
a basal tic- andesite composition) also have much less 
spectral structure in this spectral region than water 
vapor [Bandfield et al . , 2000b]. Hematite has some 
significant spectral structure at the frequencies of in- 
terest, but it is at the 50-100 cm -1 scale [Christensen 
et al . , 2000]. Hematite features are observed over only 
a very limited portion of the surface and we expect 
the features to be mostly removed by the continuum 
removal process. Therefore, the total uncertainty in 
the water vapor band depth index caused by the con- 
tinuum removal is small. 

The above four sources of uncertainty occur in the 
process of finding the water vapor band depth index. 
All four are relatively small and we estimate that the 
total uncertainty in the index for 6.25 cm -1 resolution 
data to be no more than 0.003 (with typical values 
of the index ranging from 0.03 to 0.10, see Fig. 3). 
The uncertainty in the index for 12.5 cm -1 resolution 
data may be somewhat higher because of the fewer 
number of points used in forming the index. The 
largest uncertainty in obtaining a final averaged band 
depth index for a particular L s /latitude bin is likely 
caused by uncertainty in the calibration of 12.5 cm -1 
resolution index values relative to those at 6.25 cm -1 . 
Figure 3 shows the relation between the two, and the 
spread of points around the given linear relation gives 
an indication of the uncertainty which appears to be 
roughly 0.005-0.010 (i.e. 10-20%) in a root-mean- 

square sense. 

The final two sources of uncertainty come from the 
process of calibrating the water vapor band depth in- 
dex to a column abundance in precipitable microns. 
The computation of a water vapor band depth in- 
dex based on a synthetic water vapor spectrum in- 
volves a different set of uncertainties than that de- 


scribed above for computation of the observed band 
depth. There is no longer a contribution from in- 
strument noise and calibration or from removal of 
a continuum opacity from other sources. However, 
there are new uncertainties not present when using 
the observed data, such as uncertainty in the H 2 O 
line strengths, positions, and shapes (line widths will 
be discussed separately, below) and error caused by 
use of the correlated-k approximation. The first three 
uncertainties are very difficult to estimate, but based 
on the excellent fits we get to the observed data (see 
Fig. 4), we believe this is not a dominant source of un- 
certainty. Results from the correlated-k method can 
be directly compared with “exact” line-by-line calcu- 
lations to evaluate the error involved with using the 
correlated-k approximation. With the 20 correlated- 
k integration points that we use, errors are typically 
< 1 - 2 %. 

The largest source of uncertainty in the end-to-end 
retrieval algorithm is in our assumption of a value for 
the broadening ratio, defined above as the amount 
of water vapor line- broadening by CO2 compared to 
that by terrestrial air (as tabulated in the GETS A and 
HITRAN databases). This uncertainty enters as a 
systematic error essentially multiplying all results by 
a nearly-constant factor. The value of 1.5 assumed 
here leads to an increase in the water vapor band 
depth index of about 20% over that for terrestrial 
air-broadening. The individual water vapor lines be- 
ing observed are highly saturated and are therefore in 
the regime where opacity goes as the square- root of 
abundance. Thus, abundance goes as opacity (or the 
band depth index which is based on opacity) squared 
and our nominal value for the broadening ratio gives 
an abundance ^-40% lower than would be obtained 
with terrestrial air-broadening. The likely value for 
the true broadening ratio (averaged over the impor- 
tant lines) is between 1.2 and 1.8 (see Section 3.5) 
meaning that there could be an error in retrieved wa- 
ter vapor column abundance of order ±25%. Note 
that this error would affect all retrieved water vapor 
amounts by a similar factor and would have very little 
effect on the dependence of water vapor as a function 
of L t and latitude. 

In summary, the largest source of random error is 
the cross-calibration of 6.25 and 12.5 cm -1 data. The 
total estimated random error is 0.005-0.010, which 
corresponds to an uncertainty in water vapor column 
abundance of about 3-5 pr-/im. A significant sys- 
tematic error of as much as 25% is possible through 
uncertainty in the broadening ratio. This would re- 
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suit in a nearly-constant multiplicative factor being 
applied to all water vapor column abundances and 
would not affect the distribution of water vapor as a 
function of season and latitude. 

4. Results 

The top panel of Plate 1 shows the column abun- 
dance of water vapor observed by TES as a func- 
tion of season (/**) and latitude. Water vapor abun- 
dance is shown for the MGS mapping period between 
L m = 104°, March 1, 1999 and L a = 137°, March 28, 
2001 . Results include only daytime data (local time 
~1400) where the surface temperature is >220 K. Re- 
sults have been zonally averaged into bins that are 2 
wide in latitude and 2° wide in L a . The total inte- 
grated water column integrated perpendicularly from 
the surface to infinity is shown. No attempt has been 
made here to correct for topography or to scale to an 
equivalent pressure surface. 

4.1. Overview 

The TES water vapor results (top panel of Plate 1) 
show a strong maximum in the northern hemisphere 
summer. A peak water column of ~100 pr-/im was 
attained in each of the two Martian years observed 
poleward of 80° N latitude at fs 3 — 110°-120°. At. 

= 120°, water vapor abundance dropped off mono- 
tonically from north to south with an abundance of oO 
pr-/im at 60°N latitude, 20 pr-/jm at 15°N latitude, 
and less than 5 pr-/im of water vapor south of 20° S 
latitude. After L a = 1 30° , water vapor abundance de- 
clined rapidly at high northerly latitudes failing be- 
low 10 pr-/jm by L a = 170°. However, the water 
vapor abundance at low northerly latitudes contin- 
ued to increase during this time forming a tongue 
of high water vapor abundance extending equator- 
ward and to later season from middle northerly lat- 
itudes. For any particular latitude in the northern 
hemisphere, the peak water vapor column abundance 
occurred progressively later in the season for lower 
northern latitudes. At 80°N latitude, water vapor 
column abundance reached a maximum value at L a = 
115°, as compared to L a = 135° for 45°N latitude and 
[ Jjt = 150° for 30°N latitude. The equatorial water 
vapor column exhibited a broad maximum between 
L a = 180° and 210°. By L a = 170° a well-developed 
maximum in water vapor appeared between the equa- 
tor and 30° N latitude. This maximum in water va- 
por as a function of latitude at 0°~30 o N latitudes 
persisted throughout northern fall and winter until 


f Jg — 40° in the following year, when water vapor 
abundance rapidly increased throughout the northern 
hemisphere. 

Tn the southern hemisphere there was a gradual rise 
in water vapor abundance throughout southern win- 
ter and spring as water vapor appeared to be trans- 
ported southward from the northern hemisphere. In 
late southern spring (L a = 220°) , water vapor at low 
southerly latitudes decreased somewhat, while wa- 
ter vapor at middle and high southerly latitudes in- 
creased. A maximum in southern hemisphere water 
vapor was reached at L a = 290° with a column abun- 
dance of ~40 pr-/jm poleward of 75° S latitude. At 
this time there was a local minimum in the latitudinal 
dependence of water vapor at 10 O- 30°S latitude and a 
small local maximum at 0° and 30° N latitude, which 
may be caused by transport of water vapor from the 
southern hemisphere to the north via the Hadley cir- 
culation. The latitudinal dependence of water vapor 
in the southern summer at L a — 290° was very dif- 
ferent than that in the northern summer at L $ = 
110°. While there was significant water vapor at high 
latitudes near the summer pole in both seasons, the 
northern hemisphere summer maximum was stronger 
and did not show a secondary low-latitude maximum 
in the winter hemisphere. 

Water vapor abundance decreased from the south- 
ern summer polar maximum after L$ = 300 , and 
decreased planetwide after L a — 330° ■ The period 
between L a = 330° and 40° exhibited minimum water 
vapor abundance with the average global abundance 
of 7-8 pr-/im. Water vapor abundance began to in- 
crease significantly in the northern hemisphere after 
[ Js = 40°, steadily climbing to its peak value L , = 
120° while the southern hemisphere water vapor re- 
mained relatively unchanged at a very low level. 

4.2. Comparison with Viking MAWD 

The bottom panel of Plate 1 shows water vapor 
column abundance retrieved from the Viking Or- 
biter Mars Atmospheric Water Detectors (MAWD) 
[Farmer ei al 1977; Jakosky and Farmer, 1982]. 
There are some broad similarities between the MAWD 
and TES results. The most obvious similarity is the 
magnitude and timing of the north polar summertime 
maximum in water vapor abundance. Both MAWD 
and TES show >75 pr-/im water vapor poleward of 
75° N latitude from L a = 100°-125°. In fact, the gen- 
eral dependence of water vapor on L t and latitude 
between about L a = 75° and 140° is fairly similar, 
although TES water vapor abundance is somewhat 
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higher at low latitudes. Both data sets also indicate 
that the timing of the northern hemisphere maximum 
occurs later in season with decreasing latitude. An- 
other similarity between the MAWD and TES results 
is the existence of a generally low amount of water va- 
por in middle and high latitudes in the winter hemi- 
sphere (in both the north and the south). 

Although the above features of the MAWD and 
TES water vapor L 8 /latitude maps are similar, there 
are also a number of very significant differences. Fig- 
ure 5 shows the ratio of TES to MAWD water vapor 
abundance as a function of L, and season. The largest 
distinction arises from the appearance of a strong 
maximum in water vapor at south polar latitudes just 
after summer solstice in the TES results. No such 
maximum appears in the MAWD results. Interannual 
variability likely accounts for some of this difference, 
but some of the difference is associated with an ob- 
servational bias in the MAWD data caused by two 
global, or “planet-encircling” dust storms which be- 
gan near L t = 205° and 274°. As indicated in the 
MAWD papers [ Farmer et a/., 1977; Davies , 1979; 
Jakosky and Farmer , 1982], scattering of the solar 
radiation by dust and water ice biases the MAWD 
water vapor abundance to lower values because the 
observed radiation does not sample the entire column 
of water vapor. Davies [1979] estimated that a visible 
dust opacity of 0.5 would result in a derived water va- 
por abundance 10% lower than actual. Jakosky and 
Farmer [1982] estimated that a visible dust opacity of 
1.0 would result in a derived water vapor abundance 
20% lower than actual, and they suggested that the 
MAWD water vapor abundance should be considered 
as a lower limit when visible opacity is greater than 
unity. 

During the first Viking year (1977), there were two 
planet-encircling dust storms. Drops in MAWD wa- 
ter vapor at all latitudes are apparent in Plate 1 at 
both L, = 205° and 270°, with the low values per- 
sisting beyond L t = 330°. Visible aerosol opacity 
measured by the two Viking Landers [Colburn et al. 
1989] indicates that opacities greater than unity ex- 
isted during the entire period between L t = 205° 
and 330° for both landers. The shading in Fig. 5 
shows L B /latitude bins where the infrared dust opac- 
ity inferred from the Viking Orbiter Infrared Ther- 
mal Mapper (IRTM) instrument [Martin, 1986; Mar- 
tin and Richardson , 1993] is >0.4 (equivalent to vis- 
ible opacity >~0.8). The IRTM dust opacities are 
>0.4 for the period L t — 205°-330° over much of 
the planet, in good qualitative agreement with the 


corresponding visible dust opacities of Colburn et al 
[1989]. The solid line bounding the contours on the 
top and bottom in Fig, 5 represents the latitude limit 
for which both TES and MAWD data exist (IRTM 
dust opacity coverage is comparable). As further in- 
dication of the large influence of dust on MAWD re- 
sults, the one area of clearing between the two planet- 
encircling dust storms at L t = 250°-270°, 30°-60°S 
latitude corresponds closely with a noticeable increase 
in MAWD water vapor abundance (see Plate 1). Nev- 
ertheless, it is interesting to note that although the 
MAWD results are clearly biased by dust scatter- 
ing, real interannual variability may contribute to the 
differences between TES and MAWD results. The 
Viking IRTM did not observe significantly elevated 
dust opacity at the highest southerly latitudes dur- 
ing the second planet-encircling dust storm (L s = 
280°-310\ 70°-85°S latitude) [Martin and Richard- 
son, 1993] as seen in Fig. 5. Dust opacities there 
were comparable to those in the relatively clear pe- 
riod between the two dust storms mentioned above 
(l 3 — 250°— 270°, 30°-60°S latitude). The lack of any 
increase in MAWD water vapor in the region L t = 
280°-310°, 70°-8o 0 S latitude despite the high water 
vapor observed then by TES suggests that the south- 
ern hemisphere water vapor maximum may have been 
significantly reduced in 1977. 

We conclude that for the dusty period between 
l s = 205° and 330°, 1) the MAWD water vapor re- 
sults are biased low by a as much as a factor of 3 and 
should be regarded as a lower limit of actual water 
vapor abundance, and 2) the L s and latitude depen- 
dence of the MAWD water vapor results is not ac- 
curate. This is consistent with the recommendations 
made in Jakosky and Farmer [1982]. The water vapor 
abundance retrieved from TES spectra are only very 
minimally affected by dust because dust scattering is 
much less important at TES wavelengths (^30 /im) 
than at the MAWD wavelength (1.38 /im). Further- 
more, dust loading levels during the TES operations 
(1999-2001) [Smith et al , 2001] were much smaller 
than those present during the 1977 period of MAWD 
observations. 

Outside of the dusty period L a = 205° -330° there 
remain small, but significant differences between the 
MAWD and TES water vapor maps. Most apparent 
is at low northerly latitudes (0°-30°N) between L, — 
150° and 205° where the water vapor abundance de- 
rived from TES data is roughly 20-40% higher than 
that derived from MAWD data. There is also a 
more pronounced latitudinal gradient between L $ = 
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0° and 60° in the TES results, with higher water va- 
por ahundanee at 0°-30°N latitude and lower water 
vapor ahundanee at middle southerly latitudes than 
observed by MAWD. The inerease in water vapor be- 
tween i s = 40° and 90° at 30°-70°N latitude is a 
gradual process in the TES results, while the MAWD 
results show little water vapor inerease at. these lati- 
tudes until an abrupt inerease between i s — 60° and 
75°. These differences are likely caused by interan- 
nual variability. 

Figure 5 shows that outside of the dusty i a = 
205°-330° period, the globally- averaged amount of 
water vapor is generally about 1 0—20% higher for TES 
than for MAWD. This may be an indication of in- 
terannual variability, but the difference is also well 
within the combined systematic (overall scaling) error 
bars for the TES and MAWD water vapor results. In 
the periods immediately preceding and following the 
two planet-encircling dust storms (i 9 = 180°-205° 
and i 3 = 330°-360°), the TES results are system- 
atically higher by 30-75%. This could be caused by 
interannual variability, by a small bias in the MAWD 
results during these moderately dusty times, or by a 
combination of both. 

4.3. Comparison with Ground-Based and 
Other Observations 

In addition to the Viking MAWD data, there ex- 
ists a large body of ground-based and space-based 
observations of Mars water vapor made at wave- 
lengths from the visible to microwave. The three 
main ground-based programs that have monitored 
Martian water vapor extensively over the years are 
the ~8200- A observations taken at McDonald Obser- 
vatory [ Barker , 1970; Jakosky and Barker , 1984] and 
by a group at the University of Arizona [Rizk ct ai , 
1991 ; Sprague ct a/., 1996; 2001], and the cm and mm- 
wave observations taken by Clancy ct al. [1992; 1996]. 
More recently, a series of observations has been taken 
at the Infrared Telescope Facility (I RTF) by Novak 
ct ai [2001]. Space-based observations of Martian 
water vapor have been performed at thermal infrared 
wavelengths (including some of the same lines used 
by TES) by the Infrared Space Observatory (ISO) 
spacecraft [Bergdorf ct ai , 2000] and at submillime- 
ter wavelengths by the Submillimeter Wave Astron- 
omy Satellite (SWAS) [Gurwcll ct ai, 2000]. Addi- 
tional data sources for water vapor from Mars orbit 
include thermal infrared measurements from the IRIS 
instrument on Mariner 9 [Hand ct ai, 1972], from 
the Phobos mission using solar occultations made by 


the Auguste experiment [ Bodin ct ai, 1997], from the 
Phobos infrared-imaging spectrometer (ISM) [ Rosen - 
qvist ct ai , 1992], and from the Pathfinder imaging 
system (IMP) [Titov ct ai , 1999]. 

Comparison of TES water vapor abundance with 
the above observations is difficult given the large 
differences in retrieval technique, spatial resolution, 
temporal sampling, and local time of day. However, 
some general trends do emerge when the data are 
taken as a whole. Near i 9 — 120°, nearly all ob- 
servations consistently show a large peak in water va- 
por abundance at high northerly latitudes with lower 
amounts of water at low latitudes and very little wa- 
ter in the southern hemisphere. This general behav- 
ior is consistent with TES and MAWD results. The 
L t range of 340°-90° (northern spring) seems to be 
a period with a relatively large amount of interan- 
nual variability. Results of Ftxzk ct ai [1991], Clancy 
ct ai [1992], and Sprague ct ai [1996] vary by as 
much as a factor of three during this period. TES 
water vapor abundance tends to fall near the up- 
per limit of this range of variation for the northern 
hemisphere, and near the middle or lower limit for 
the southern hemisphere. Near-infrared ground-based 
data from the IRTF taken by Novak ct ai [2001] at 
i t = 67° (January 1997) is consistent with TES data 
taken two Mars years later southward of 30°N lat- 
itude, but is significantly lower at middle and high 
northerly latitudes. Data from the Phobos/TSM in- 
strument [Rosenqvist ct ai, 1992] taken during this 
season (9 pr-/im for 0°~20°S latitude, and ~12 pr- 
pm for 0°-20°N latitude at i s = 2°-18*) are in good 
agreement with TES data. 

At i s = 150°-160°, both the Pathfinder [Titov ct 
ai, 1999] and the SWAS [Gurwcll ct ai, 2000] find 
very low water vapor abundances of 6 and 8 pr-/jm, 
respectively. TES observed 24 pr-pm water vapor at 
the time and latitude of the Pathfinder observation 
(is = 150°, 19°N latitude), and a disk-integrated 
value of about 20 pr-pm at the time of the SWAS 
observation (i s = 160°). It is not clear why these 
two observations are so low compared to TES (and to 
MAWD). 

There are relatively few additional observations 
that can be used to determine the degree to which 
the southern summer maximum observed by TES at 
i 9 = 260° -320° repeats each Martian year. Barker 
[1970] and Jakosky and Barker [1984] present disk- 
integrated water vapor from McDonald Observatory 
taken during this season. Following Jakosky and 
Barker [1984], we estimate a disk- in teg rated TES 
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value of 18 pr-pm water vapor at L s — 220° increas- 
ing to 23 pr-pm at L $ — 280°-300° and falling to 16 
pr-/im by L s = 330°. The McDonald Observatory 
data from 1969, often cited as an unusually wet year 
are roughly consistent with the TES values through 
L, — 300°. However, the very high values of 30- 
45 pr-/im observed at L a = 320°-350° are at least 
twice as high as those observed by TES at that sea- 
son. The extensive McDonald Observatory measure- 
ments made during 1972-1973 are consistently some- 
what higher than MAWD, but agree better with the 
TES data, especially for L, — 150°-300°. Hand et al. 
[1972] report 10-20 pr-/im water vapor from Mariner 
9 IRIS spectra at L t = 297° and 32 1 0 in the region 
around the south polar cap, and observations by Rizk 
et al. [1992] made in December 1988 at L s = 320° 
show a water vapor abundance of 1 5-20 pr-/jm at high 
southerly latitudes. These values are lower than those 
observed by TES (30-40 pr-pm at L t = 297°and 20- 
25 pr-/im at L t — 320°), but substantially higher than 
the 5 pr-/im observed by MAWD. 

The best ground-based measurements for compar- 
ison are those taken since March 1, 1999 so that a 
direct comparison can be made with the TES water 
vapor retrievals. We know of three groups that have 
observations of Martian water vapor during this time. 
The program at McDonald Observatory to monitor 
Mars water vapor led by E. Barker has observations 
from 1999 and 2000 [personal communication, 2001] 
but these data have not yet been reduced and cali- 
brated. Sprague et al. [2001] report on observations 
taken at 70°-80°N latitude between at L, = 105°- 
115° and at 150° in 1999 which can be directly com- 
pared with TES observations taken from the begin- 
ning of mapping. TES water vapor abundance at 70°- 
80°N latitude and L, — 1 05° —1 1 5° cover the range 
~50-100 pr-/jm and so are roughly consistent, though 
perhaps somewhat higher than the 43-76 pr-pm re- 
ported by Sprague et al. [2001]. At 70°-80°N latitude 
and L s = 150°, TES observations give 20-30 pr-/jm 
water vapor which is significantly higher than the 9- 
17 pr-/jm reported by Sprague et al. [2001]. Prelim- 
inary results are available from recent observations 
by R. Novak, M. Mumma, and M. DiSanti [personal 
communication, 2001] taken in the near-infrared (3.6 
/im) at the TRTF on March 22, 1999 (L, = 113°) 
and January 16, 2001 (h s = 104°). The 1RTF water 
vapor results are comparable to those of TES from 
the same date and latitude near the equator, but are 
significantly lower (by a factor of 2 or 3) than TES 
at middle and high northerly latitudes. This is the 


season where TES (and others) have consistently ob- 
served the northern hemisphere summertime maxi- 
mum of water vapor. 

5. Discussion 

5.1. Condensation Level 

An important quantity that comes out of the re- 
trieval process is the water condensation level. The 
condensation level is computed in the retrieval algo- 
rithm to set the vertical distribution of water vapor. 
Figure 6 shows the water condensation level in terms 
of kilometers above the surface as a function of L s and 
latitude. The condensation level is a function of atmo- 
spheric temperature and water vapor abundance, but 
the dependence on temperature is the dominant con- 
trolling factor. The latitude and seasonal dependence 
in Fig. 6 largely follows that of averaged atmospheric 
temperatures [compare with the top panel of Plate 1 
of Smith et ai } 2001]. Between L s = 0° and 180° 
the atmosphere is relatively cool and the condensa- 
tion height is low (10-20 km above the surface). As 
atmospheric temperatures increase toward perihelion 
(currently at L t = 251°) due to enhanced solar flux 
and dust loading [Smith et a/., 2001], the condensa- 
tion level rises to altitudes above 40 km. The abun- 
dance of water vapor plays a relatively minor role, 
with the condensation level inversely related to water 
vapor abundance. This is most clearly seen at places 
of very low water vapor, such as 40 C> -50°S latitude 
from L t = 60°-210° where the condensation level in- 
creases to the south despite warmer temperatures to 
the north. 

The inferred seasonal variation of the water con- 
densation level is large. Water vapor is essentially 
well mixed throughout the atmosphere for half the 
year ( L, — 1 80°— 360°) , but confined to the lowest 
scale height or two for the other half of the year (L e =r 
0°-180°). This large, seasonal variation in the con- 
densation level is just one more component of the two 
distinct seasons that appear to exist on present-day 
Mars, with relatively warm and dusty conditions in 
the perihelion (southern spring and summer) season, 
and relatively cool, dust-free, and cloudy conditions 
in the aphelion (northern spring and summer) sea- 
son. The seasonal behavior described above for the 
condensation level was emphasized by Clancy et al. 
[1996] who suggested that the seasonal variation in 
the height of the water condensation level could play 
a pivotal role in the transport of water vapor and the 
observed asymmetry in the latitudinal dependence of 
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water vapor near solstice. Near perihelion the con- 
densation level is relatively high and water vapor can 
be transported from south to north via the Hadley cir- 
culation without condensing. Near aphelion the con- 
densation level is relatively low and water vapor can 
not be transported via the Hadley circulation with- 
out condensing. This mechanism could explain the 
relative maximum observed at low northerly latitudes 
during the fall and winter, and the much larger latitu- 
dinal gradient in water vapor in the northern summer 
as compared to the southern summer. 

5.2. Globally- Averaged Trends 

The top panel of Fig. 7 shows the total amount of 
atmospheric water vapor as a function of season in- 
tegrated over the northern and southern hemispheres 
and over the entire globe. Condensed water is not 
included in the total integrated amount of water be- 
cause its contribution is insignificant. To perform 
the integration over an entire hemisphere, an esti- 
mate was made for water vapor abundance in the 
high-latitude regions where no results are available 
because of insufficient thermal contrast between the 
surface and the atmosphere. In the fall and win- 
ter seasons in each hemisphere (L s = 180°-360° in 
the northern hemisphere and 0°-180° in the southern 
hemisphere), a low nominal value of 3 pr-/im water 
vapor was used to be consistent with the low val- 
ues retrieved at the highest latitudes where good re- 
sults were obtained during those seasons (see Plate 
1). In the spring and summer seasons, a constant 
water vapor abundance was extrapolated to the pole 
from the most poleward retrieved results. These ex- 
trapolations have little effect on the hemispherically- 
and globally- aver aged water vapor amounts for two 
reasons. Firstly, the surface area where water vapor 
is extrapolated (which is near the polo) is a small 
fraction of the total surface area over a hemisphere. 
Secondly, the water vapor values being extrapolated 
are generally small (except during early summer) and 
contribute little to the total. 

The total global amount of water vapor varies from 
~1 .0-2.5xl0 15 g (equivalent to 14-35 pr-/im or 1 .1-2.7 
km 3 ice) during the year. The major seasonal feature 
in the global amount of water vapor is a ~50% drop 
from L b — 310° -360° and a corresponding recovery 
from L s = 50°-120°. Both the northern and south- 
ern hemispheric totals show maxima in the spring and 
summer seasons and minima in the fall and winter. 
The annually- aver aged total amount of water vapor 
is 1.2xl0 15 , 0.7xl0 15 , and 1.9xl0 15 g (equivalent to 


17, 9.5, and 13 pr-/im) for the northern hemisphere, 
southern hemisphere, and the entire globe, respec- 
tively. 

As expected, the total water vapor amounts shown 
in Fig. 7 are slightly higher from L t = 330° -205° than 
those computed by Jakosky and Farmer [1982] and by 
Habcrlc and Jakosky [1990] using MAWD results, and 
significantly higher during the dusty L t = 205°-330° 
season. Unlike the nearly constant level observed by 
TES, MAWD results show a significant decrease in 
the global amount of water vapor between L , = 170° 
and 320° because of the two planet-encircling dust 
storms that occurred during 1977 at that time. 

The changes in total atmospheric water vapor as a 
function of season give clues to the important large- 
scale sources, sinks, and transport processes for wa- 
ter vapor. Between L s = 50° and 120°, northern 
hemisphere water vapor increases dramatically while 
southern hemisphere water vapor remains at a nearly 
constant low level. The increase in northern hemi- 
sphere water vapor occurs primarily at middle and 
high latitudes (see Plates 1 and 2) and so the likely 
sources for the water vapor are the receding seasonal 
ice cap and the release of water adsorbed into the 
regolith. The north residual polar cap can also con- 
tribute to the increase in northern hemisphere water 
toward the end of this period after = 80° when 
it is exposed. The relative importance of the seasonal 
and residual ice caps versus the regolith as a source 
(or as a sink) can not be determined from this analysis 
alone. 

The northern hemisphere peak in water vapor 
abundance lasts for only a short time. During the 
period L s = 140°-240°, the amount of northern hemi- 
sphere water vapor decreases rapidly to a value less 
than half of its peak value at L B = 120°. Although 
global water vapor shows a small (~20%) decline from 
J Jt = 140°-170°, it is essentially constant from L t = 
170°-320° because of an increase in southern hemi- 
sphere water vapor which largely offsets the decrease 
in northern hemisphere water vapor. This suggests 
that both non-atmospheric sinks and transport con- 
tribute to the loss of northern hemisphere water va- 
por. During the first portion of this period (before 
Lg = 180°), the loss is dominated by incorporation of 
water back into the residual and seasonal ice caps and 
by adsorption of water into the regolith. Later, there 
also appears to be a net transport of water vapor from 
the northern to the southern hemisphere that is most 
active between L t = 150° and 230°. 

The amount of northern hemisphere water vapor 
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stops declining near L a = 230°, but the amount of 
southern hemisphere water vapor continues to in- 
crease somewhat until reaching a peak at L a = 290°. 
During this period, water vapor equatorward of 40° S 
latitude slowly decreases while water vapor poleward 
of 40°S latitude rapidly increases. The latitudes 
where water vapor increases are clear of seasonal CO 2 
frost by this time of the year [Christensen et ai y 
2001], so the source for this increase seems to be pri- 
marily from the release of water adsorbed into the 
regolith and possibly from transport. Beginning at 
L $ = 310°, southern hemisphere water vapor declines 
rapidly without a compensating increase in the north- 
ern hemisphere water vapor, such that the northern 
spring season presents a significant minimum in the 
annual variation of Mars global atmospheric water. 
This contrast in global water abundance between the 
northern and southern spring seasons is a major as- 
pect of the north-south asymmetry in atmospheric 
water behavior, as described in section 5.3. There 
appears to be little, if any, transport from the south- 
ern hemisphere to the northern hemisphere during the 
decay of the southern summer water vapor maximum. 
The decline in southern hemisphere water vapor ap- 
pears to be dominated by losses to the seasonal cap, 
water vapor being cold-trapped onto the residual CO 2 
polar cap, and adsorption into the regolith. 

5.3. Hemispheric Asymmetry 

The bottom panel of Fig. 7 shows the total amounts 
of water vapor in the northern and southern hemi- 
spheres as a function of season in terms of A L 9 after 
spring equinox so that the behavior of the two hemi- 
spheres can be directly compared. There arc signif- 
icant asymmetries. The overall level of water vapor 
throughout the year is higher in the northern hemi- 
sphere. The amount of northern hemisphere water 
vapor never reaches the very low levels that south- 
ern hemisphere water vapor does in the fall because a 
moderate amount of water vapor (at least 10 pr-/im) 
is maintained at low northerly latitudes during the en- 
tire year (see Plate 1 and further discussion in section 
5.4). Both hemispheres exhibit very similarly- timed 
early-summer maxima in water vapor and then large 
decreases in water vapor in late summer and early fall. 
However, the amplitude of the northern hemisphere 
maximum is larger, and the time-dependence of the 
increase in water vapor leading up to the summer 
maxima are very different for the two hemispheres. 
In the north, the increase begins rather abruptly at 
about A L a ~ 40° and the entire increase occurs be- 


tween A L a = 40°and 120°. In the south, there is 
a similarly-timed increase in water vapor (A L a = 
40°-120°) but it is much smaller. The majority of 
the increase in southern hemisphere water vapor oc- 
curs earlier during the winter and early spring and 
appears to be caused by transport of water vapor 
from the northern to the southern hemisphere. Be- 
cause of these differences, in the early northern sum- 
mer [L a = 100°-140°), the northern hemisphere holds 
~85% of the total global atmospheric water, while 
in the early southern summer (L s — 280°-320°) the 
southern hemisphere holds only ~60% of the total 
global atmospheric water. 

Significant asymmetry also currently exists in the 
seasonal and latitudinal distribution of water va- 
por between the northern and southern hemispheres. 
This asymmetry is shown in Plate 2. The water va- 
por abundances shown in Plate 2 are the same as 
those shown in the top panel of Plate 1 except that 
the southern hemisphere data have been “flipped” in 
latitude (so that south pole is at the top and the equa- 
tor is at the bottom) and shifted by 180° of L a rel- 
ative to the northern hemisphere data so that a di- 
rect comparison between the two hemispheres can be 
made. There are several key differences. As men- 
tioned above, the summertime maximum that is ap- 
parent at high latitudes in both hemispheres is about 
twice as strong in terms of maximum water vapor col- 
umn abundance in the north as it is in the south. 
A second key difference is the existence of a deep 
minimum in water vapor covering the entire southern 
hemisphere during the fall and early winter. No such 
minimum exists in the northern hemisphere. In fact, 
there is more water vapor at low northerly latitudes 
during the fall and winter than there is in the spring. 
The third and perhaps most striking difference be- 
tween the two hemispheres is the completely different 
behavior of the decay of the summertime maximum 
between mid-summer and mid-fall. Tn the northern 
hemisphere, an area of relatively high water vapor 
extends increasingly equatorward as the season pro- 
gresses. Near northern fall equinox the water vapor 
appears to cross the equator and continue southward 
eventually contributing to the formation of the sum- 
mertime maximum in water vapor at high southerly 
latitudes. In stark contrast, the southern hemisphere 
maximum rather abruptly decays in middle to late 
summer with no indication at all of an area of rela- 
tively high water vapor extending equatorward, and 
little indication of northward transport of water va- 
por. 
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5.4. Latitude and Longitude Dependence of 
Water Vapor 

In all of the above analysis we have always aver- 
aged over longitude to view the dependence of water 
vapor with L, and latitude. However, retrieved water 
vapor column abundance values can also be binned in 
latitude and longitude (averaging over time) to pro- 
duce time-averaged latitude-longitude maps. The top 
panel of Plate 3 shows such a map where the time 
average is taken over the entire annual cycle. Im- 
mediately apparent is an inverse correlation with to- 
pography. Hollas shows relatively large water vapor 
columns while Tharsis shows relatively low water va- 
por. This kind of inverse correlation implies that the 
water vapor mixing ratio is constant or a relatively 
weak function of altitude in the lowest scale height or 
two where most of the mass of the atmosphere lies. 
Thus the retrieved results are consistent with our al- 
gorithm assumption that water vapor is well-mixed 
up to the condensation level which is typically 1-1.5 
scale heights for L t = 20°-150° and higher for the 
rest of the year. 

The effect of topography can be removed by divid- 
ing retrieved water vapor abundance by a quantity 
proportional to surface pressure. The bottom panel 
of Plate 3 shows annually- aver aged water vapor di- 
vided by the quantity (pnurf/51), where p 8ur f is the 
annually-averaged surface pressure in mbar. Surface 
pressures are computed using the topographic map 
produced by the MGS Mars Orbiter Laser Altimeter 
(MOLA) team [Smith at at. , 1999] along with the an- 
nual variation of surface pressure [Hess et at., 1980; 
Tillman ct ai, 1993]. A map of annually- aver aged 
surface pressure is shown in Fig. 8. We will use the 
term “scaled water vapor” for the average annual wa- 
ter vapor abundance divided by (p» u rf/01)- With this 
correction, the dependence of water vapor with lati- 
tude and longitude that is independent of topography 
is now apparent. 

Figure 9 shows the latitudinal dependence of the 
scaled water vapor obtained by taking a zonal aver- 
age of the results shown in Plate 3. The hemispheric 
asymmetry at middle and high latitudes is caused by 
hemispheric asymmetry in altitude. When the lower 
northern hemisphere altitudes are taken account of, 
the northern and southern hemispheres poleward of 
45° latitude are nearly identical. This is a markedly 
different result than that obtained by Jakosky and 
Farmer [1982] using the MAWD data who found a 
large gradient in annually-averaged water vapor with 
twice as much water vapor at high northerly latitudes 


than at high southerly latitudes (even after correcting 
for the hemispheric topographic difference). 

The most obvious feature in the latitude depen- 
dence of scaled water vapor (Fig. 9) is a pronounced 
maximum between 10°S and 40°N with roughly 50% 
more water vapor in this latitude band than outside it. 
The annually-averaged water vapor abundance is high 
in this latitude band because water vapor is main- 
tained year-round at these latitudes (see Plate 1). 
Water vapor may be maintained at low northerly lat- 
itudes during late northern fall and winter by trans- 
port from the southern hemisphere via the Hadley 
circulation. The northern boundary of the band of 
high water vapor (~40°N latitude) is likely defined 
by the southernmost extent of the north seasonal po- 
lar cap. Equatorward of that latitude the sinks for 
water vapor are much reduced and water vapor can 
remain at a moderate level year-round. 

The bottom panel of Plate 3 shows that there are 
large longitude variations in annually-averaged scaled 
water vapor equatorward of ~50° N and S latitude. 
The maximum seen in Fig. 9 (and Plate 1) at 10 S- 
40°N latitude is primarily caused by two large re- 
gions of high water vapor abundance from 60°-l 50°W 
longitude (Tharsis) and 300°-10°W longitude (Ara- 
bia Terra). Although the highlands around Olym- 
pus Mons and the Tharsis volcanoes have a high 
scaled water vapor, the volcanoes themselves do not. 
Large longitude variations exist outside of the lati- 
tude band of high annually-averaged water vapor as 
well. At 25°S latitude, scaled water vapor varies from 
its global minimum of 7.5 pr-/im at 35° W longitude 
(just south of Margaritifer Sinus) to a value twice as 
high just 90° of longitude to both the cast and west. 

The top panel of Fig. 10 shows the albedo of 
Mars derived from TES visible bolometric observa- 
tions [Christensen et of, 2001]. The low latitude 
longitude dependence in scaled water vapor shows 
a strong positive correlation with albedo (high wa- 
ter vapor in high albedo regions). This same cor- 
relation with albedo was also noted by Jakosky and 
Farmer [1982] using MAWD results and by Rosen- 
qvist et al. [1992] using results from the ISM instru- 
ment on the Phobos spacecraft. At low latitudes wa- 
ter vapor is also anti-correlated with surface pressure 
[Smith et ai , 1999] (Fig. 8), and with thermal in- 
ertia [Mellon ct ai , 2000] (bottom panel of Fig. 10). 
All three of these quantities (albedo, surface pressure, 
and thermal inertia) have moderate correlations (or 
anti-correlations) between them and so it is difficult 
to sort out the cause and effect. Albedo, thermal 
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inertia, and surface pressure (topography) have di- 
rect relations to the average and diurnal variation of 
surface and subsurface temperatures, the pore space 
size and grain size of the regolith, and to the cir- 
culation patterns of the atmosphere. These can in 
turn influence the transport of water vapor and its ad- 
sorption and desorption into the regolith. It is likely 
that a combination of albedo, thermal inertia, and 
topography contribute to the longitude dependence 
of annually- averaged water vapor. Poleward of 40°N 
latitude the correlation of scaled water vapor with 
albedo, and the anti-correlation with thermal inertia 
and surface pressure, abruptly vanishes, and perhaps 
even changes sign. In the southern hemisphere, the 
correlation of scaled water vapor with albedo persists 
to at least 60°S latitude. Poleward of that there is 
little longitude dependence in the scaled water vapor. 

Scaled water vapor might still be anti -correlated 
with surface pressure even after the topography has 
been “removed” because water vapor is not com- 
pletely well-mixed through the entire atmosphere. 
Therefore, dividing by (p sur f/6.1) can actually be a 
slight overcompensation for topography. Support- 
ing this idea is the fact that the correlation between 
scaled water vapor and surface pressure is lower in 
the southern hemisphere than in the northern hemi- 
sphere. For example, Hellas and Argyre do not show 
up at all in the bottom panel of Plate 3 (scaled wa- 
ter vapor) although Hellas is very prominent in the 
top panel of Plate 3 (unsealed). In the southern 
spring and summer when southern hemisphere water 
vapor abundance is high (and the most contribution 
to the annual average is made), the condensation level 
is high and water vapor is more nearly well-mixed 
throughout the atmosphere. In the northern spring 
and summer when northern hemisphere water vapor 
is high the condensation level is lower and water vapor 
is not as completely well-mixed. 

6. Summary 

We have developed an algorithm to retrieve water 
vapor column abundance from TES nadir-geometry 
spectra using rotational lines from 240- 360 cm -1 . 
Observations of water vapor by TES during the Mars 
Global Surveyor mapping mission now cover over one 
full Martian year and allow an analysis of the de- 
pendence of water vapor on latitude, longitude, and 
season (L s ) in unprecedented detail. 

The TES results show a maximum in water va- 
por abundance at high latitudes during mid-summer 


in both hemispheres. The southern hemisphere max- 
imum was not observed by Viking MAWD. Results 
from MAWD taken between L s - 205° and 330° ap- 
pear to be greatly affected by two planet-encircling 
dust storms. Scattering by dust and water ice biases 
the MAWD water vapors to lower values since the ob- 
served radiation does not sample the entire column of 
water vapor. As stated by Jakosky and Farmer [1 982] , 
we believe that MAWD water vapor results from this 
period must be regarded as lower limits. The TES- 
derived water vapor results are only very minimally 
affected by dust scattering because scattering is much 
less important at TES wavelengths (~30 ft m) than 
at the MAWD wavelength (1.38 /Jtm), and no planet- 
encircling dust storm was observed during the TES 
year (1999-2000). Outside of the dusty L s = 205°- 
330°season, TES and MAWD are generally compara- 
ble, although TES results are systematically 10-50% 
higher than MAWD results, and significant differences 
are apparent in the details of the dependence on L 9 
and latitude. 

There are large differences in the hemispheric (north 
versus south) and seasonal (perihelion versus aphe- 
lion) behavior of water vapor. The summertime peak 
in water vapor abundance is roughly twice as high in 
the north as in the south. The decay of the northern 
hemispheres summertime maximum leads to a rela- 
tive maximum in water vapor that extends further 
southward with time, eventually crossing the equator 
to the southern hemisphere. The decay of the south- 
ern hemisphere summertime maximum has no corre- 
sponding northward transport. Water vapor abun- 
dance declines monotonically at all latitudes at the 
end of southern summer. 

Primarily because of the seasonal dependence of 
atmospheric temperatures, there is a strong seasonal 
dependence on the water condensation level. The 
condensation level is typically 10-15 km (1-1.5 scale 
heights) above the surface during the aphelion season 
(k a = 0°-180°), but greater than 30 km (^3 scale 
heights) during the perihelion season (/,, = 180° 
360°). This means that water vapor is confined to 
lower altitudes in the aphelion season, and transport 
between the two hemispheres seems to be influenced 
by this. A very steep latitudinal gradient in water 
vapor abundance forms during the aphelion season 
(high in the north), while water vapor is spread more 
uniformly in latitude during the perihelion season. 

The annually- averaged amount of water vapor is 
17 pr-/im in the northern hemisphere and 9.5 pr-pm 
in the southern hemisphere. However, when the ef- 
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feet of topography is removed, the annually- averaged 
amount of water vapor becomes 17 pr-/im (refer- 
enced to a 6.1 mbar pressure surface) in the latitude 
band from 10°S-40°N, and 12 pr-/im everywhere else. 
Northerly and southerly high latitudes have the same 
scaled water vapor abundance. This latitudinal de- 
pendence is very different than the analysis of M AWD 
results by Jakosky and Farmer [1982] which showed 
a large gradient from north to south in annually- 
averaged water vapor abundance even after correction 
for topography. The latitude and longitude depen- 
dence of annually-averaged water vapor retrieved by 
TES (corrected for topography) has a significant pos- 
itive correlation with albedo and significant negative 
correlations with thermal inertia and surface pressure. 

Additional analysis and observations promise to 
make the TES data set an even more useful and com- 
prehensive tool for the study of the Martian water cy- 
cle. Continued observations by TES during the Mars 
Global Surveyor extended mission will provide the 
amount of interannual variability between two con- 
secutive Martian years, and further analysis of the 
existing TES spectra taken in the limb-geometry will 
provide crucial information about the vertical struc- 
ture and diurnal (~0200 versus ~1400 local time) de- 
pendence of water vapor. 
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Figure 1 . A synthetic line-by-line spectrum of water vapor convolved to the TES resolution of 6.25 cm -1 for the 
wavelengths used in this analysis for typical conditions (solid line). The spectrum is shown in terms of opacity 
The filled boxes show the best-fit wavenumber locations for the TES samples in the 6.25 cm" 1 resolution mode. 
The numbers refer to TES sample numbers. The five water vapor band peaks used in this analysis (centered near 
TES samples 20, 26, 30, 35, and 40) are separated by roughly 5 TES samples or 25 cm -1 . 

Figure 2. A comparison of averaged TES spectra at 6.25 cm"' (solid line, filled boxes) and 12.5 cm" 1 (dotted 
line, open boxes) spectral resolution taken at the same season and latitude (L, = 306°-308°and 78°-80°S latitude). 
Spectra arc shown in terms of opacity. Water vapor bands are much better defined at 6.25 cm-’. 

Figure 3. The relationship between water vapor band depth computed using TES data taken in the 12.5 cm -1 
resolution mode (t , 2 5 ) and that taken in the 6.25 cm -1 resolution mode (r 6 . 25 ). A nearly linear relation exists 
allowing band depth at 12.5 cm -1 resolution to be converted to an equivalent band depth at 6.25 cm -1 resolution. 
The linear fit is given by r s 25 = 1.45 (r 12 5 - 0.0035). 

Figure 4. An example of the fit between the observed averaged TES spectrum (open boxes and solid line) and the 
best-fit synthetic spectrum (stars and dotted line) for one particular /,,/latitude bin (centered at 43°N latitude, 
L, - 85°) shown in terms of opacity. An opacity continuum has been removed from the observed spectrum. A 
water vapor column abundance of 25 pr-/jm is indicated. 

Figure 5. Contour map of the ratio of water vapor abundance derived from TES spectra to the water vapor 
abundance derived from MAWD data [Jakosky and Farmer , 1982] as a function of L s and latitude. The shad- 
ing represents L t and latitudes where Viking IRTM measurements showed 9-/im dust opacity >0.4 [Martin and 
Richardson, 1993]. The solid lines bounding the contours on the top and bottom represent the latitude extent of 
where both TES and MAWD results exist. Planet-encircling dust storms occurred at — 205° and 270°. The 
times and locations of high dust opacity (shaded areas) correspond to the times and locations of high TES/MAWD 
water vapor ratio. 

Figure 6. The height in km above the surface of the water condensation level as a function of L s and latitude. 
Contour levels higher than 40 km above the surface are not shown. The solid lines bounding the contours on the 
top and bottom represent the latitude extent of where results exist. 

Figure 7. The top panel shows total amount of atmospheric water vapor as a function of season for the northern 
hemisphere (dashed line), the southern hemisphere (dotted line), and the entire globe (solid line). The bottom panel 
shows the hemispheric data in the top panel replotted to allow direct seasonal comparison of the two hemispheres. 
Season is expressed as AL 3 from spring equinox ( L s = 0° for the northern hemisphere and = 180° for the 
southern hemisphere). Equivalent values are 1 km 3 of ice is 0.917xl0 15 g, and 10 pr-um over a hemisphere is 
0.723x1 0 15 g. 

Figure 8. Map of the annual average value of surface pressure in mbar based on MOLA data [Smith at a/., 1999]. 
In this map, quantities have been smoothed by the same amount as the contours shown in Plate 3. 

Figure 9. The latitudinal behavior of the annually-averaged water vapor abundance, and the latitudinal behavior 
of the annually-averaged water vapor abundance divided by (jw/6.1) to remove the effect of topography. 

Figure 10. Maps of (top) broadband visual albedo, and (bottom) thermal inertia in J/m 2 K a 1 / 2 . In these maps, 
quantities have been smoothed by the same amount as the contours shown in Plate 3. The albedo [Christensen e.t 
al, 2001] and the thermal inertia [Mellon ct ai, 2001] maps are from TES measurements. 
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Plate 1. The column abundance of water vapor as a function of I., and latitude: (top) as observed by TES^ 
Contours show a smoothed representation of the results, and (bottom) as observed by Viking M AWD [Jakosky and 
Farmer, 1982]. 


Plate 2. The distribution of water vapor as a function of season and latitude for the (top) northern and (bottom) 
southern hemispheres. Season is expressed as A/„ from spring equinox (L, = 0° for the northern hemisphere an 
L, = 180° for the southern hemisphere). Both hemispheres are shown with the pole at the top and the equator a 
the bottom for direct comparison. Contours show a smoothed representation of the results. 


Plate 3. Maps of (top) seasonally-averaged water vapor column abundance and (bottom) seasonally-averaged water 
vapor column abundance divided by (jw/6.1) to remove the effect of topography. Contours show a smoothed 
representation of the maps. 
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